AD*A117 828 NAVAL POSTGRADUATE SCHOOL MONTEREY CA F/G 9/2 

AN INTERACTIVE COMPUTER PROGRAM FOR THE PRELIMINARY DESIGN AND ==ETC(U) 
. MAR B82 J t PAQUETTE. 

UNCLASSIFIED 


NAVAL POSTGRADUATE SCHOOL 


Monterey, California 


THESIS 


AN INTERACTIVE COMPUTER PROGRAM 
FOR THE PRELIMINARY DESIGN AND 
ANALYSIS OF MARINE REDUCTION GEARS 
by 
Joseph Louis Paquette 


March, 1982 


Thesis Advisor: G. Cantin 


Approved for public release; distribution unlimited. 


Seaeteneterepecsncromanenaminbormanan on n.ano_socniniaaaaanamm annie manmecnmtsinr mance eigenen aaa rane aT 2D ena lapaatateaatemted 


UNCLASSIFIED 


SECUMTY CLASSIFICATION OF THIS PAGE (When Date Entered) 


REPORT DOCUMENTATION PAGE 


2. GOVT ACCESSION NO. 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


5. TYPE OF REPORT & PEMOO COVERED 
Master's Thesis 
March 1982 


6. PERFORMING ORG. REPORT NUMBER 
. Con GR GRA 


4. TITLE (and Subsite) 
An Interactive Computer Program 


for the Preliminary Design and 


Analvsis of Marine Red on Geg 
+ AUTHOR a) 


Joseph Louis Paquette 


wT. PROJECT Tasn 
tT NUMBERS 


12. REPORT CATE 
March 1982 
13. wmumsGER OF PAGES 
187 


18. SECURITY CLASS. (of thie répert) 


Naval Postgraduate School 
Monterey, California 93940 


. CONTROLLING OF FICE NAME AND ADORESS 


Naval Postgraduate School 
Monterey, California 93940 


UNCLASSIFIED 


a. OECL ASSIFICATION/ DOWNGRADING 
SCwEQULE 


Approved for public release; distribution unlimited. 


17. CISTRIGUTION STATEMENT (of the abatroct entered in Bleck 20, !f different trem Report) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue an reverse side if nececcary and identity by bieck number) 
Computer Aided Design; Marine Gear Design; Marine Gears; Helical 


Gears; Computer Aided Marine Gear Design; Conventional Parallel 
Axis Gears; Epicyclic Gears; REGAD; Computer Code 


20. AGSTAACT (Continue en reverse ofde if neeecoary end identify by bieck mambor) 


av 
The objective of this project was to develop an interactive 
computer program providing flexibility in the design and analysis 


of marine propulsion gears. The program, Reduction Gear 
Analysis and Design (REGAD), will handle conventional parallel 
axis and simple epicyclic reduction gears. It is capable 9 


DO ,7on'7y 1473 carrion oF 1 wov e613 OBsoL ETE 
S/™ 0102-084* 6601 | SECUMTY CLASSIFICATION OF THIS PAGE / we Entered) 


UNCLASSIFIED 


‘ eer EE A LEE ET REE ERE SL TE TI IE EES I 
| Racumiy Ci 2660016 211ON OF Y0Is POGE/Men Nose Bnsornd 


generating preliminary designs of new gear sets or of providing 
analyses of existing or proposed gear sets. Program development, 


organization, and operation are discussed. . 
ry 


Pistributions ~ 


AVALlability Codes | 


Dice Rtas and/or 
“"  " Speciaz 


| DD Form, 1473 UNCLASSIFIED 
ot s/n 0102-014-6601 2 SUCUATY Ci, AEmMPICATION EF Tutt PaGEGren Bare Basoved: 


|. 12 erg REO RERUN 


Approved for public release; distribution unligited. 


An Interactive Computer Prograa 
for the Preliminary Design and 
Analysis of Marine Reduction Gears 


by 


Joseph Louis Paquette 
Lieutenant, United States Navy 
B.S.S.E., Onited States Naval Academy, 1976 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
March, 1982 


Authors 


Thesis Advisor 


YMA 


Second Reader 


Chairman, Depar 


Dean of Science and Engineering 


3 


ABSTRACT 


The objective of this project was to develop an interac- 
tive computer program providing flexibility in the design 
and analysis of marine propulsion gears. The prograa, 
Reduction Gear Analysis and Design (REGAD), will handie con- 
wential parallel axis and simple apicyciic reduction gears. 
It is capable of generating preliminary designs of new gear 
sets or providing analyses of existing or proposed gear 
sats. Program development, organization, and operation are 


discussed. 
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Ie JZNTRODUCTION 


In the conceptual stage of ship design, many parameters 
and options are considered. This is especially true with 
respect to the propulsion plant. Changes in hull design, 
displacement, and numbers of propellers all affect changes 


in the requirements for the propulsion plant. There are 


also various options under consideration in the propulsion 


plant: turbines or internal combustion engines, the number 
of engines, the auxiliaries required for support, etc. All 
of these will affect che initial design of the reduction 
gears. It is, therefore, important to be able to produce 
reliainary designs of reduction gears for the options under 
consideration. 

Preliminary designs provide useful information on feasi- 
ble arrangements and size without going into the specific 
design details dependent upon manufacturing. Since any nun- 
ber of preliminary designs may be required due to perturba- 
tions discussed above, it is important to be able +o 
automate the design process. An interactive computer pro- 
gram providing various options would free the engineer from 
*edious, «ime consuming, and often error prone number amanip- 
ulation and allow him to produce multiple designs for con- 


sideration. I* would also provide a quick means of checking 


she effects of various parameters in addition to the ability 
+o analyze proposed designs or configurations. 

Reduction Gear Analysis and Design (RZGAD) was developed 
to fill this need. It is an interactive computer program 
offering close user control through numerous options. Being 


interactive, it provides a rapid means of designing or ana- 


lyzing a gear set, thereby reducing the turn-around time 


inheren* in the use of batch systems. The program was kept 
modularized and well documented for ease of maintenance and 
modification. The modularized construction also provides an 
additional banefit of being able to use this program on 


smaller computers by using an overlay scheme. 


II. PROGRAS CAPABILITIES 


A. SCOPE 

REGAD was written to provide preliminary designs or 
analyses of marine propulsion reduction gears. It is inca- 
pable of providing detailed designs or performing detailed 
analyses since specifics of manufacturing are not reguired 
for input. The program does not consider shafting, bear- 
ings, lubrication, couplings, casings, or other auxiliaries. 
It will provide sizing information in the form of pitch 
diameters, 2ffective facewidths, gear ratios, and numbers of 
teeth per gear. In addition, the program will provide esti- 
mates of loadings and stress levels. Estimated weight and 
dimensions of the gear set are also provided. 

All computations are based primarily on the American 
Gear Manufacturers Association's standards [Ref. 1, 2, 3] 
using appropriate constants for marine propulsion gears 
(Ref. 4, 5]. As an option in the program, these constants 
can be replaced by the user to enable him to investigate 
other applicaticns such as reduction gears for ships service 


or emergency generators. 
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B. LIMITATIONS AND OPTIONS 
Program application is limited to marine reduction gears 
with a maximum of three reduction stayes. Conventional par- 
allel axis and simple epicyclic arrangements with helical 
gears are possible. when dealing with epicyclics, it is 
assum2d that load sharing of the planets is achieved and 
that the ring gear is suitably flexible. Efficiencies of 
the gear sets are not provided since power losses are not 
computed. While estimates of bending and contact stresses 
are provided, scoring can not be estimated since lubrication 
is not considered. REGAD does not require the K-factor as 
input as in previous programs since hardness ranges for the 
pinions and gears are required. However, the K~factors are 
computed and displayed for reference purposes. The weight 
estimates are based on actual designs and do rot include 
*urning gears, attached lubrication oil pumps, or other 
auxiliaries. 
The following is a list of major options provided by «he 

program: 

(1) brief, on-line program description 

(2) choice of design or analysis 

(3) listing of preprogrammed constants and an ability to 

change selected constants of the user's choice 
(4) choice of single, double, cr triple reductions 
(5) choice of single or double helical gears 


(6) choice of six hardness ranges for gears and pinions 
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(7) conventional parallel axis arrangements (see Figures 
1, 2, and 3) 
(a) one or two power inputs 
(b) single power path (articulated) or dual power 
paths (locked train) 


(8) simple epicyclic arrangements (see Figure 4) 


| 
| 
| 


(a) choice of planetary or star arrangements 
(b) single power input 


(c) choice of three, four, or five planet/star gears. 
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IIT. PROGRAM ORGANIZATION AND OPERATION 


A. REGAD FLOWPATHS 

As stated previously, the program was designed in nodu- 
lar form with aach module consisting of a number of subpro- 
grams. These aodules are just conceptual groupings of 
associated subprograms, and are not related to actual pro- 
gram implementation on any specific computer. Figure 5 
shows the basic flow paths of the program. Module One is 
for program initialization and problem set up. Module Two 
performs calculations for conventional parallel axis gear 
sets, while Module Three handles epicyclic gear sets. Mod- 
ul2 Pour is a grouping of all the computational subprograms 


required by the other modules. 


B. MODULE IDENTIFICATION AND DESCRIPTION 
This section provides a brief description of each sub- 
program in 2ach module. 
1. Module One : Initialization and Set-up 
Module One contains the subprograms necessary for 
initialization, execution, and initial data entry. It is, 
basically, the control module for the program. The follow- 


ing is a grouping of the subprograms in Module One. 
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REGAD 
REGAD is the main program. It provides the 
options for either design or analysis and either parallel 
axis or epicyclic arrangements and controls the flow to tae 
proper module. I* then calls the required subprograms for 
execution of Modules One, and Two or Three. 
b. BLOCK DATA 
The BLOCK DATA subprogram initializes variables 
in each of the common blocks. 
c. SUBROUTINE DSCRPT 
This subroutine is called by REGAD after an 
affirmative response to a user option to provide a brief 


description of the REGAD package. It contains an option to 


stop «he program if only a program description is desired. 


d. SUBROUTINE INPUT 
All options and initial design parameters 
entered via “his subroutine which is called by REGAD. 
2. SUBROUTINE AGMA 
The constants for marine propulsion gears 
required by various AGMA formulations are initialized 
BLOCK DATA subprogram, and can be listed as an option in 
REGAD. REGAD calls this subroutine after an affirmative 
user response +o display the preprogrammed values. This 
subroutine then allows the user to selectively change any 


constant desired. 


2. Module Two : Parallel axis 

This module contains all the major subprogranas 
called by REGAD to provide an initial design or to perform 
an analysis of convential parallel axis reduction gears. 
The following is a grouping of the subprograms in Module 
Two. 

a. SUBROUTINE PRLDES 

This subroutine will produce a design of a par- 

allel axis gear set. All pinion and gear diameters, effec- 
tive facewidths, and gear ratios are computed using a basic 
random search optimization technique to find a feasible 
design by attempting to minimize a function of gear pair 
volume. I+ should be noted that, while attempting to mini- 
mize gear volume, the design is not necessarily optiaized 
for minimum weight. The optimization technique is used here 
only to produce a feasible design in terms of dimension and 
power constraints by minimizing a function of gear pair vol- 
ume. To produce a truly optimized design for mininua 
weight, a full optimization must include many more design 
variables such as helix and pressure angles, pitches, and 
hardnesses in addition to the dimensions. Additional con- 
straints such as stress and unit load levels would need to 
be incorporat2d. All of this would require a more sophisti- 
catad and efficient optimization technique than is used 


here. 
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b. SUBROUTINE PRLANL 
To analyze a proposed or existing design, REGAD 
will call this subroutine. It will request, as user-sup- 
plied input, the basic information calculated in PRLDES, 
i.e., pitch diameters and effective facewidths. Using this 
information, PRLANL will compute other parameters such as 
gear ratios, power and speed splits, and numbers of teeth 
per gear. 
c. SUBROUTINE PRLRES 
Immediately following a call to PRLDES or 
PRLANL, REGAD will call PRLRES to compute all remaining 
information such as expected loadings and stress levels. 
The user should be aware that the stress levels are computed 
accoriing to AGMA formulations {Ref. 2 and 3] and take into 
accoun* load distribution and overloads. This will produce 
levels *ha* may seem high but are actually closer to actual 
levels to be expected in service. 
ad. SUBROUTINE PRLSIZ 
REGAD calls this subroutine after PRLRES *o con- 
pute estimates of gear set weight and gearbox dimensions. 
These estimates are determined by empirical relationships 


obtained from a rather limited data base of actual designs. 
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e. SUBROUTINE PRLOOT 

This is the last subroutine called by REGAD in 
the parallal axis path. It provides a detailed output of 
the results obtained from the design or analysis including 
design parameters entered by the user, the dimensions of 
2ach component, expected loadings and stress levels, and 
configuration information. 

3. Module Three : Epicyclic 


Module Three contains all the major subprograags 


called by REGAD to design or analyze simple epicyclic reduc- 


tion gears. Subroutines EPCDES, EPCANL, EPCRES, EPCSIZ, and 
2PCOUT are all analogous to those in Module Two. They pez- 
form «he same functions, but for simple epicyclic gears. 
Therefore, individual descriptions will not be repeated 
here. 
4. Module Four : Computational Subprogram Library 
This module is an organizational grouping of all the 
subprograms called by those in Modules One, Two, and Three. 
a. Subroutine Subprograas 
The following are the subroutines used: 
(1) GFI - subroutine to compute the AGMA dur- 
ability geometry factor, I 
GFJ - subroutine to compute the AGMA 


strength geometry factor, Jd. 


Real Punction Subprograrzs 


The following are the real function subprograms 


(1) ARCCOS - computes the arc cosine of two 


arguments 


ARCSIN ~- computes the arc sine of two argu- 


nents 


AGMAE1 - uses LaGrangian interpolation of 


Table E-1 (Ref. 1] to compute the constants 


required for the stress concentration factor 


foraulation 


(4) CKDATA - called by SUBROUTINE AGMA to allow 


the user to change the preprogrammed con- 


stants 


(5) POWERB - computes allowabie service power 


based on AGMA strength rating (Ref. 3] 


(6) POWERH - computes allowable service power 


based on AGMA durability rating [Ref. 2] 


RTFNDR - a modified version of FUNCTION 


ZEROIN (Ref. 6] used to find a zero of a 


function in a specified interval 


FALFA - the function reguired by SUBROUTINE 


GFJ and the zero af which is computed in 


FUNCTION RTFNDR 


SHRLD - computes the load sharing ratio, ny 


(10) THICK - computes tooth thickness at any 
diameter given a known thickness at a dif- 


ferent diameter. 


C. DATA TRANSFER 
All data transfer between subprograms in Modules One, 

Two, and Three is via combinations of seven common blocks. 
Data transfer to and from subprograms in Module Four is via 
argument lists and common blocks as required. The following 
is a list of the common blocks used: 

(1) /AGMAB/ : constants for AGMA strength formulations 

(2) /AGMAH/ ; constants for AGMA durability formulations 

(3) /DESDAT/ : design parameters and options 

(4) /SDESPRL/Y : parallel axis design information 

(5) /SRESPRL/ :; parallel axis computational results 

(6) /DESEPC/ : epicyclic design information 

(7) /RESEPC/ : epicyclic computational results. 
The variables in each common block along with their defini- 


tions can be found in Appendix B. 


D. PROGRAM OPERATION 

REGAD is an interactive program designed to allow the 
user to solve his problem at a terminal. Being interactive, 
*he program has many options that control program execution, 


in addition to requests for data necessary for the execution 
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of the program, Each request for information will contain 
the necessary guidelines needed by the the user to respond. 
This may take the form of a mini-table containing informa- 
tion on each option choice, the range of values when a spe- 
cific quantity is requested, or units, where applicable, of 
+he raquested data. 

All option parameters are integer values and should not 
be antered with a decimal. Option codes entered by the user 
are checked for validity to ensure they fall within the 
allowed range. If two options are offered, enter a 1 ora 
2. <Any value entered less than one will automatically 
default to one, and any value greater than two will automat- 
ically default to two. In those cases where there are more 
than «wo options, the response is checked to see if it falls 
within the allowed range. If it does not, a message alerts 
the user to this fact and allows him to re-enter the correct 
code. Some questions require affimative or negative respon- 
ses. To reply, use a Y for yes or an N for no. Use of 
cher values may give undesirable results. 

Every attempt has been made to anticipate possible error 
conditions. If one of these is encountered, a message is 
generated to inform the user. If tha error encountered is a 
*erminal error, the message will also indicate that the pro- 


gram tun was aborted under program control. 
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| A detailed development of this package is provided in 
| Appendix A where specifics can be found. Appendix 3 pro- 
vides a cross-reference of the variables used in Appendix A 
with «hose used in the program. It also contains detailed 
information on the common blocks. Sample runs of the pro- 
gram can be found in Appendix C, and a complete listing of 


the program is in Appendix D. 
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IV. CONCLUSIONS AND RECO N 

Computer aided design (CAD) is an important and useful 
tool for engineers. AS computer technology continues to 
expand, CAD will become increasingly available for the prac- 
sicing engineer, allowing him to use his initiative in 
design instead of being a slave to the numbers involved. 
REGAD is such a tool for use in the preliminary design of 
marin? reduction gears during the conceptual stages of pro- 
pulsion plant design. 

REGAD could become even more useful if additional 
options ate provided. A module to perform sensitivity anal- 


yses of a given design would greatly enhance the use of this 


program. This option would allow the user to start with any 


dasign and vary a selected variable over a specified range 
to determine its impact on the design. It could also he 
used to "fine tune" a design by modifying selected parame- 
ters to produce the results desired without having to rerun 
the orogram for each nodification. Graphics would add 
another dimension by providing graphical displays of the 
gear arrangements and of certain data such as the results of 
a sensitivity analysis. A module to handle various compos- 
ite designs of parallel axis and epicyclic gears would be an 


importan+ addition. Also, it is recommended that a larger 


data base be collected to provide more accurate empirical 


constants for the weight and gearbox size estimates. 


APPENDIX A 


PROGRAM D 


With the exception of several qeneral conversion rela- 
tionships, all computations are accomplished in Modules Two 
and Three with calls +o subprograms in Module Four. The 
analytical relationships used in the program will be exan- 
ined, however, most of the relations used can be easily 
found in the literature and in various texts, so background 


developments will not be given. 


I. GeNERAL RELATIONS Hes 


The following relationships are used in Module One and 
in various other subprograms. The transverse diametral 
pitch of any gear is the ratio of its number cf teeth to its 


pitch diameter; 


N 
,=— (1) 
d 


The nornal and «ransverse diametral pitches are related by; 


P,= Picos ¥ (2) 


and <=h2 pressure angles by; 


tan ®, = tan @, cosy (3) 


Axial pitch is defined as; 


——— (4) 


II. CONVENTION ARALLEL AXIS FORMULATION 
Subroutines PRLANL and PRLDES each provide the pitch 
diame+ters of the pinions and gears, the effective face- 


widths, the stage reduction ratios, the numbers of teeth per 


gear, speed and power splits, and the geometry factors to 
subroutines PRLRES and PRLSIZ to compute all further infor- 
mation. Thea speed splits are the actual speeds of the indi- 
vidual gears and a power split is the actual power 
transferred by a gear. The strength and durakility geometry 
factors are computed in separate subroutines in Module Four 
\ 


and will be discussed later. 


Power splits are determined from the configuration. For 
a single power path configuration, the power is transferred 
aqually from «he pinicn to the gear, where in a dual power 


A. COMMON RELATIONSHIPS 
| pata configuration, the pinion transfers one half its power 
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to each of two gears. These splits are computed exactly 
since losses are neglected. 
Speed splits and stage reduction ratios are based on; 
D Ng 
Bez? = (5) 
d Dp 
Numbers of teeth on each gear are computed from the 


equation below and are rounded to the nearest integer. 


N =dx PR, (6) 


B. DETERMINATION OF DIAMETERS AND FACEWIDTHS 

All diameters and facewidtns are entered by the user in 
subroutine PRLANL. Stage gear ratios, power and speed 
splits and aumbers of teeth per gear are computed as dis- 
cussei in the previous section. 

In subroutine PRLDES, the diameters, facewidths, and 
stag2 gear ratios are determined by using a basic local ran- 
dom search optimization technique to produce a feasible 
design. This algorithm requires an initial design to start. 

The initial design is based on Dudley's [Ref. 7] formu- 
lation for preliminary estimates of gear size; 


31500 Pwr  (mg*1)3 
C2PF = 


(7) 


K Np Dg 
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a 
c === (Me +1) (8) 
2 


By substituting equation 8 into equation 7, a formula for 
estimating pinion diameter is obtained; 


126000 Pwr (Ng + 1) 


——— (9) 


np K (F/d) Rg 


where F/d = 1.0 for single helical gears and F/d = 2.25 for 
doubl2 helical gears. The term K is the K-factor which is 
an indication of durability. An expression for estimating K 
is provided by Thoma [Ref. 4]; 


Ss. x 10-*\2 /2.80 


Ce oes oF 


(10) 


wher2 =the constants used are tne AGMA durability constants. 
The K-factor in equation 10 is for the second reduction. 
Por the first reduction, multiply K from equation 10 by 


1.20, The initial estimates for the stage gear ratios are: 
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(1) single reduction a,= M 
(2) double reduction mg=VM,+ 3 dual power path 


Bg =V4,- 1 single power path 


i Mg, = VE, /Bo, 
| (3) triple reduction ag,= VE, 
! mg,2 VS, + 3 
a= Vi /2g2q, - 
The initial facewiiths used are: 
(1) single helical gears F=d 
(2) double helical gears PF = 2.25 4d. 
With this initial design as a starting point for «he 
random search algorithm, successive designs are determined 


by randomly adding small amounts of between 1.0 and -1.0 *o 
*he diameters, facewidths, and stage gear ratios. These 

| small amounts are scaled to take into account the difference 
in ctange of values for each variable. This process will 
attempt to find a feasible dasign in which all specified 
constraints are satisfied. If the initial design violates 
one or more constraints, the design that violates them the 
least in succeeding iterations will be kept until a design 
satisfying ali constraints is found. Once a feasible design 


is found, an attampt to improve this design is made by try- 


\ 
q 
! 

ing to reduc2 the size of the gears by minimizing a function 

of gear pair volume; 
| 
Voiume = DQ cer = YY Cilmget)2a2P] (11) { 
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0 Re. 


The interior summation is over the number of reduction 


stages, and =he exterior summation is over the nuaber of 


power inputs. The constraints imposed which determine the 


limits on each of tne designs are: 


(1) actual transmitted power is less than or equal to the 


allowabla service power in accordance with references 


2 and 3 


(2) maximum gear diameter of 200 inches due to manufactur- 


ing limitations 


(3) minimum facewid=h greater than four axial pitches to 


ensure proper helical action 


(4) maximum facewidth less than the pinion pitch diameter 


for a single helical gear or 2.25 times the pinion 


| pitch diameter for a double helical gear 


(5) pinions and gears in succeeding reduction stages are 


to be larger than those in the previous stage due to 


the greater amounts of torque carried 


(6) ia dual power path arrangements, the gear ratio for 


stage due to the torgue carried. 


The design obtained can than be adjusted by the user as 


desired by changing parameters with the analysis option. 


i 
| 
4 
| 
| 
each sraduction stage is greater than the preceeding 


Tit. C ¢€ RMULATION 


As in Modul? Two, the pitch diameters, effective face- 
widths, stage reduction ratios, numbers of teeth, speed and 


power splits, and the geometry factors are all entered or 


computed in the analysis or design subroutines (EPCANL or 

EPCDES) for use in the final computations subroutines 

(EPCRES and EPCSIZ). Here, the speed splits are the rota- 

tional speeds of the sun and planet gears and of either the 

ting gear or the carrier, depending on the configuration. 

Planetary arrangements have fixed ring gears while star 

arrangements have fixed carriers. Also, the direction of 

| rotation must be considered. Star arrangements reverse the 

| direction of rotation of the input and the planetary 

arrangements will maintain direction of ro-ation. Assuming 


equal load sharing of the planets and neglecting losses, 


power splits are straightforward. The input and output pow- 


ers are equal while each planet carries an equal share of 
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the total power. Load sharing is an important consideration 
in the design of epicyclic gears, and must be assured in 
marin2 reduction gears due to the high power levels experi- 
enced. Equal load sharing of the planets can be reasonably 


achieved in several different ways. One method requires the 


sun gear +o float, supported only by che planet gears, with 


a relatively flexible ring gear to allow for inaccuracies in 
the teeth. There are also mechanical devices available to 
assist in achieving an equal division of the load. Experi- 
ence has shown, for marine applications, that three to five 
planets with stage ratios in the range of two to eight work 


best. 


A. COMMON RELATIONSHIPS 

Unlike conventional parallel axis arrangements, there 
ar2 specific numerical rules governing the proper assembly 
an@ operation of an epicyclic gear set. These involve the 
selection of «he numbers of teeth and planets along with 
computing the various speed ratios. Mesh frequencies are 
also configuration dependent as seen in a following section. 

There ar2 basically three relationships that must be 
satisfied to ansure proper assembly and operation. The 
firse is a r2lationship defining the speed ratio of the epi- 
cyclic stage since it is not nerely the ratio of numbers of 
seeth or diameters as in a conventional parallel axis gear 
sex. The second relationship requires the ring gear diame- 
ter to be aqual to the sum of the sun gear diameter and 
swice a planez gear's diameter. This ensures the planets’ 
ability to Fit between the sun and ring gear. For the final 
Telationship, it can be shown geometrically that the sum of 


the numbers of *eeth on tne sun gear and ring gear must be 
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an integral multiple of the number of planets in the gear 
set to ensure proper alignment and meshing of all teeth. I[t 
should be noted tha+ these relationships are based on 
equally spaced planets around the sun gear. The above rela- 
+ionships are conveniently expressed in terms of numbers of 
teeth on gach gaar as seen in references 7 and 8. The speed 


ratio for a planetary arrangement is; 


Ro Np 
n, N, 


and for a star arrangement; 


=n, -Ny 
io =e: (13) 
nj Ng 


wher? the negative sign indicates the stir arrangement's 
teversal of rotational direction of the input. The rota- 


sional speed of the planet gears is required for the design 


of their bearings and can be determined by; 


Ren ao (14) 


Noun 
where no in each equation above is the speed of the carrier 
for a planetaty arrangement or is the ring gear's speed for 
a star arrangement. The assesbply and meshing relations in 


*“erms of tooth numbers are; 
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Na = No + 2 Now (15) 
ands; 
Na, + Ne = k NP (16) 


where k is an integer and NP is tne number of planets. 


B. DETERMINATION OF DIAMETERS AND FACEWIDTHS 

Th2 random search technique discussed for conventional 
parallel axis gears is used to provide the epicyclic diame- 
ters, facewidths, and stage gear ratios. Equation 9 is used 
so provide an initial estimat2 of sun gear diameters where m, 
is replaced by <he ratio of the planet's pitch diameter to 


the sun's pitch diameter. This value is usually in the 
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range of 1.5 *o 3; therefore, a random number in this range 
is used to stars the problem. Once the sun gear diameter is 
estimated, «he other diameters can be found using the rela- 
cionships in equations 12 to 16. The initial estimates for 
the stage gear ratio are the roots of «he overall ratio cor- 
responding *o the number of reduction stages. For example, 
Mg, aad Rg, for a double reduction gear set would be the 
square root of the overall ratic. Initial facewidths are 
chosen as before. The initial estimates of the diameters, 


facewidths, and the gear ratios provide a starting point for 
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the random search algorithm discussed previously. Again, 
the method will attempt to improve feasible designs by sini- 


mizing a function of gear volume; 


Volume = 5° (NP-d2 + d2 + d2)-F (17) 


where the summation is over the number of reduction stages. 

The constraints imposed are similar to those for the paral- 

lel axis gears: 

(1) actual transmitted power is less than or equal to the 

allowable service power in accordance with references 
2 and 3 
aaximum ring gear diameter of 150 inches due *o manu- 
facturing limitations 
minimum facewidth greater than four axial pitches to 
ansure proper helical action 
maximum facewidth less than the sun's pitch diameter 
sor a single helical gear or 2.25 times the sun's 
pitch diameter for a double helical gear 
planet gears are to be larger than sun gears due to 
she gr2ater amounts of torgue carried 
stage gear ratics are to be between 2 and 8 for each 
reduction stage. 

As before, once a design is obtained, the user can utilize 


the analysis option to obtain the desired results. 
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Once the geometry is determined in-the analysis or 
design subroutines, the couputational results subroutine 
(PRLRES or EPCRES) and the size estimates subroutine (PRLSIZ 
or EPCSIZ) are called to provide design information concern- 
ing tooth loads, stresses, and other configuration, geon- 
atry, and size information. This section describes the 
formulations used. 

The facewidth to diameter ratio is computed using the 
effective facewidth and the pitch diameter of the pinion for 
parallel axis gears or the sun gear for epicyclics. Center 
distance is taken as the average of the pinion and gear 
pitch diameters. A center distance is computed for epicycl- 
ics by finding the average of the sua and a planet gears! 
pitch diameters. 

Pitchline velocity, V, is determined by; 

wd Np 
(18) 
12 
wher2 V is in fee« per minute, dis in inches, and n, is in 


revolutions per minute. The tangential component of tooth 


load, W,, is computed from; 


126000 Pwr 
4, = (19) 
np d 


here W, is in pounds-force, Pwr is in horsepower, and d and 


Rp, are as before. Tooth loading per inch of facewidth is 


computed from; 


Tooth Load per Inch = W, /F {20) 


and the unit lead, a normalized value of the load per inch 


above, is; 


Unit Load = (21) 


where the unit load is expressed in pounds~force per square 
inch. 

Mesh frequencies rfrovide information on how often a 
tooth is loaded. Mesh frequencies for parallel axis gears 
are determined by; 

(22) 
60 
with £ expressed in Hertz. For epicyclic gears, the follow- 


ing are used: 


(a) {&°e=-___"_———" Ng (a) fs = NP Ng 
Nat Ny 
Na Ns Ng 
Neun Nat Ns Nein 
NPN, Ng 
Nat Ny Ne 


whare (a) through (c) are for planetary arrangements and (d) 
+hrough (f) are for star arrangements. 
The K-factor is computed for reference purposes by; 


Ww, (@,+1) 


‘oee,C_—-—~-—""_"- (24) 


Fd Rg 


The con*ac* stresses are computed according to reference 2 


bY; 


(25) 


_ > (26) 


Individual torgues, T, are found by; 
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d 


2000 


while the total output torgue is computed by; 


63 SHP 
(28) 
Dp 
where T has the units of thousands of inch-pounds-force ina 
both cases. Shaft horsepower, SHP, is the total power 
transferred to the oucput shaft. 
@eight and size estimates are based on empirical rela- 
*ions obtained from a limited number of actual designs. The 


relations used 


C1-( >> (d2F) Je2 
c3:})° F 


c4&-D 


Height c5:-D 


where «he constants used are found in Table 1. All dimen- 
Sions are in inches and the weight is in pounds~-force 


rounded to three significant figures. 


Table 1: Empirical Constants for Weight and 
Size Formulations 
1196.0 0.905 
0.34 0.89 
2.26 2.85 
1.20 1.30 
1.37 -- 
1, 28 1.20 


Bull Gear Ring Gear 
Diameter Diameter 


first C4: for single power inputs 
second C4: for double power inputs 


V. COMPULTATIONA UBPROGRA R ORMULATIO 
The formulations provided below are for the major compu- 
tational subprograms in Module Four. Those that are self- 
axplanatory or are not computational in nature are only 


described in general. 


A. ARCCOS AND ARCSIN 


These function subprogrags find the arc cosine and are 


sine, zespec*ively, for any two arguments. They were added 


for convenience since not all compilers have them as iater- 


nal functions. 


B. AGMAE1 

This function subprogram returns the value of the con- 
stants H, L, and M required for the determination of the 
stress concentration factor, K,, according to reference 1, 
for use in computing the strength geometry factor, J. Tapble 
E-1 in reference 1 provides the tabulated data necessary to 
perform a LaGrangian interpolation for each constant for a 
specified normal pressure angle in degrees. The interpola- 


tion formula used is; 


(O,-20) (0, -14.5) (Q,~ 14.5) (D,~25) 
F(Q,) = as 2 
57.75 -27.50 
(9,714.5) (@,-20) 
+ F, (30) 


52.50 


where F represents «he appropriate values of H, L, or 4M. 


C. CKDATA 

FUNCTION CKDATA is called by subroutine AGMA to allow 
the user to selectively change the preprogrammed constants 
by checking if the value entered is zero. If it is zero, 
she current value of the specified constant is not changed. 
This provides for flexibility in changing constants with 
multiple values, and it guards against inadvertantly enter- 


ing a value of zero. 
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D. POWERB AND POWERHA 

These function subprograms are used to compute the maxi- 
mum allowed service power, in horsepower, that can be trans- 
mitted by a gear according to references 2 and 3. The 
formulation based on the strength rating is; 


n 4d MK, F J See Ky 


Pp So eee ee 


126000 SF Kz Km Ke Py Ka Ky 


(31) 


and the durability rating formulation is; 


n a I Cy Sy, da Cu. Cy | 4 
126000 SF cy, C, Co Cm Cp Ch Cy 


wher? J and I are the respective geometry factors, F is the 
effective facewidth, n is the speed of d in revolutions per 
minut2, and d is the pinion pitch diameter for parallel axis 
or is the sun pitch diameter for epicyclics. All other val- 


ues are the oreprogrammed constants. 


E. RTFNDR AND FALPA 

The function subprogram RTFNDR, a slightiy modified ver- 
sion of FUNCTION ZEROIN (Ref. 6], is used to find the value 
of the root of the equation programmed in function PALFA. 


This root is raquired by the subroutine GFJ for the computa- 


tion of the strength geometry factor, J. 


F. SHRLD 
This function subprogram computes the load sharing ratio 
used in computing the geometry factors. The load sharing 


ratio, m,, is determined by; 


Pw T cos 0, 


(33) 
95 2 -95 2 Pag 


where Z is the length of actior defined as; 


The subscripts on the pitch diameters are: 
(1) 0 : outside diameter; d, = d + (2/P,) 
(2) b : base diameter; d, = d cos 9, 
For epicyclics, replace the outside diameters is equation 34 


with inside diameters : 4d, = d - (2/P,). 


G. THICK 
FUNCTION THICK returns the value of the normal arc 
thickness of a *ooth at a specified diameter given a *hick- 


ness at another diameter. For external gears; 
t, = d,((t,/d,) + inv 0,- inv ®, ) (35) 
and for internal gears; 


*, = 4, ((t,/d,) - inv, ¢ inv @,) (36) 


where the subscript 2 represents the desired point ard 
subscript 1 represents the known point. The involute func- 


| 
| tion is defined as; 

inv x = tan x - x (37) 
\ 


The arguments of the involute functiors in equations 35 and 
36 ar2 the transverse pressure angles at the points under 
consideration. The pressure angle at the desired point is 
defined as; 


d; cos Q, 
cos Q2= 


(38) 
d, 


The known point is usually taken at the pitch circle where 


d, = 4d, 9,=9,, and t, is defined as 


cos (39) 


ct 
w 
u 


H. GFI 
This subroutine is used to compute the AGMA durability 


geometry factor, I, in accordance with reference 2. The 


geometry factor is defined as; 


cos Q, sin 0, rN 


I = — (40) 
zc my (m, +1) 


where m, is computed by function SHRLD described above. 
plus sign applies +o external gears and the minus sign 


applies to internal gears. 


I. GFd 
SUBROUTINE GFJ is used to compute the AGMA strength 
geometry factor, J, in accordance with reference 1 with one 
major difference: the values used are from analytical devel- 
opments and are not scaled +o a normal diametral pitch of 
one as are the values used in a graphical layout discussed 
in reference 1. The strength geometry factor is defined as; 
Y~ cos? y 
(41) 
Ky Ry 
The load sharing ratio, m,, is compuczed in FUNCTION 
SHRLD. The stress concentration factor, K,, is determined 


fron; 


(42) 


H, L, and ® are determined in FUNCTION AGMAE1. The 


of the root fillet radius, r,, is; 


(b - ry)? 
(d/2cos? ~) + (b - r,) 


with the dedendua, b = 1.25/P,, and the root tip radius, 
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ry 20.28/P.. j$%The values of t and h are determined from the 
anlytical geometry of the tooth form layout described below. 
The tooth form factor, Y , is defined as; 


cos Q,, 1.5 tan Din re 


cos Q, x Cc, t 


7] 


Yo = Pra 


(44) 


where + and x are also from the tooth form layout mentioned 
previously. The helical factor, C,, is defined as; 
v 1 


i“ (45) 
100 100 


where tan v = tany sin@®, for ws 50°. The normal load 
pressure angle at the tip of the tooth, 9,,, can be seen in 


figures 6 and 7 and is given by; 


a, to 
0,,= cosmt[(———} + —— (46) 
do d, 


where the subscript o pertains to the point on the outside 
diameter and subscript b pertains to the base circle. The 
plus sign applies to internal gears and the ainus, to exter- 
nal. The thickness, to, at the outside diameter is deter- 
mined by function THICK. For internal gears, replace the 
outside values with the inside values as before. 

The graphical tooth form layout is a method by which the 


variables h, t, and x can be determined from actual 


a7 


4 
) 


measurements of a tooth form drawn and scaled for a normal 
diametral pitch of one for the case where tooth loading is 
at the tip. Loading at the tip of the tooth is the general 
practice for considering loads on helical gears. Refer to 
Figure 7 for the meanings of h, x, and T where t = 2T. 
Before determining h, x, and T analytically, several refer- 
ence parameters must be determined as suggested by McIntire 
and Lyon (Ref. 9]. The first is the radius from the center 
ef the gear to the tip of the inscribed Lewis stress para- 
bola which is pcint E in Figure 7. This point is the inter- 
section of the lin2 of action of the tip load, tangent to 
the base circle, with the tooth centerline. The radius to 


“his point is; 


d, 


dy 
T=z™ = (47) 
2 


2cos Din 

An additional reference point is required to fix the geon- 
etry. The center of of the root fillet is taken as this 
point which can be obtained by a very close approximation. 
To locate this point, the gear center is taken as the origin 
of a cartesian coordinate system with the tooth centerline 
as the vertical axis. Two possible cases exist for the 
location of this point with respect to the base circle. 
Figure 8 shows *he case where the point is inside the base 


circl2 andi Pigure 9 skows the case where it is outside. The 
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coordinates of this peint, (XC,YC), can be found from Fig- 
ures 3 and 9. For both cases it can be seen in Figures 8 


and 9 that; 


HYP = da ¢ 2, (48) 
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Figure 6, the angle, 


te 
€ = inv@ *¢ sin-! —— (49) 
a 


where t. is the chordal tooth thickness given by; 


t3 cos2 y 
ee ed Oe ea (50) 


and t+ is the nermal are tooth thickness defined earlier. 


For th2 case in Figure 8; 


XX = (HYP) sin ¢€ 


xc = XX + ©, (a) 
(51) 
Varp2- xc2 (b) 


and for the case in Pigure 9; 


re 
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(d,/2) 


@,= cos~t 


t 

| 

t 

| HYP 

{ 
: OPP, = (HYP) sin Q, 
| 
| OPP, = OPP, - 12, 
t 
| | HYP, =VoPP2 + (d,/2) 2 
! (d,/2) 
Q,= cos ~! 


A=@, t inv®, - 9, 


"=" for external gears 
"#" for internal gears (see Figure 11) 


6= At €E 
XC = (HYP) sin 6 (a) 
(52) 
YC = {HYP) cos 6 (dD) 


With *he reference values of r,, XC, and YC determined, 
she values of h, t=2T, and x can be analytically determined. 


From Pigur2 7; 


XT = ©, cosa@ (53) 
YY = :, sing (54) 
h =r, ~ YC + YH =r, - YC + ©, sina (55) 
T = (*/2) = XC - XT = XC - zc, cosa (56) 


| 
! 
{ 
| 
| 
| 


T 


YK = (57) 
tan @ 
where q@ must be determined such that; 
YK = 2h or YK - 2h = 0 (58) 
Substituting equations 53 through 57 into 58 yields; 
F(a) = XC - r,cosq@ =- 2tanq@ (r, -YC + r,sina) = 0 (59) 


Equation 59 is *he function in FALFA called by RTFNDR to 
solve for a. Once a is determined, h can be determined from 


equation 55 and T and t from equation 56. To obtain x, 


observe the following; 


y= tan-! (h/T) 
Y= (772) - Y 
and 
x = T tan, (60) 


While not precise, the identical methodology is used for 
internal gears. Figures 10 and 11 apply. The expressions 


for internal gears are given without further development; 
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-r, + YC + xr, sing 
xC ~ ©,cos a 

2T 

tan~! (h/T) 

(™/2) -¥ 


T tan Y, 


The values for h, t, and x are now usad to determine the 


stress concentration factor, equation 42, and the tooth fora 


factor, equation 44, required to compute the strength geon- 


atry factor, J, in equation 41. 


APPENDIX 3B 


LIST OF PARAMETERS 


While it is not practical to list all variables used in 
she formulations or the program, it is useful to provide a 
list of the major variables with a cross-reference between 
the analytical names and the FORTRAN names. A detailed 
listing of gach common block is also useful when studying 


+he program. 


T. DARAMETER CROSS-REPERENCE 


This section provides a listing of parameters with both 


their analytical and FORTRAN names. 


Math FORTRAN Variable 
Symbol Name Definition 

K, AKL life factor 

Km AKM load distribution factor 
Ky AKO overload factor 

Ka AKR reliability factor 

K, AKS size factor 

Kr AKT temperature factor 

Ky AKYV dynamic factor 

SF SP8 service factor 

c CDE center distance (theoretical) (in) 
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CF 

CH 

CL 

CM 

co 

CP 

CR 

cs 

ct 

cv 
SFH 
DHELIX 
HELIX 
DPHI 
PHT 
DPHIN 
PHIN 
DG 


DP 


(E=epicyclic, P=parallel axis) 
surface finish factor 
hardness factor 

life factor 

load distribution factor 
overload factor 

elastic properties factor 
reliability factor 

size factor 

*empoerature factor 

dynamic factor 

service factor 

helix angle (deg) 

helix angle (rad) 

transverse pressure angle (deg) 
transverse pressure angle (rad) 
normal pressure angle (deg) 
normal pressure angle (rad) 
diameter of gear (in) 

diameter of pinion (in) 
diameter of planet gears (in) 
diameter of ring gear (in) 
root diameter of a gear (ina) 


diameter of sun gear (in) 


facewidth (in) 


(E=epicyclic, P=parallel axis) 
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fd ratio (facewidth/diameter) 

(E=epicyclic, P=parallel axis) 

duribility geometry factor (pinion) 
GI durability geometry factor (sun) 


GEOMJG strength geometry factor (gear) 


GEOMJP strength geometry actor (pinion) 


GJS strength geometry factor (sun) 
GJPL strength geometry factor (planet) 
KFCTRE computed k-factor 

KFCTRP (E=epicyclic, P=parallel axis) 

MFE mesh frequency (Hz) 

MFP (E=epicyclic, P=parallel axis) 
overall reducticn ratio 
(Z=epicyclic, P=parallel axis) 
stage reduction ratio 
(E=epicyclic, P=parallel axis) 
number of teeth, gear 
number of teeth, pinion 

NPLNT number of planet gears in epicyclic set 

NPLN number of teeth, planet 

NR number teeth, 

NS number of teeth, 

PD transverse diametral pitch 
normal diametral pitch 
pitch line velocity (fpm) 
(E=epicyclic, P=parallel axis) 
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Pwr PWRE power split per gear pair (hp) 
PWRP (E=epicyclic, P=parallel axis) 
Din RPMIN source speed input (rpm) 
| Dout RPMOUT output shaft/propeller speed (rpm) 
nj RPMI stage input speed, epicyclic (rpm) 
| No RPMO ‘stage output speed, epicyclic (rpm) 
: Rew RPMPL planet speed, epicyclic (rpm) 
| Dp, Ng RPMP stage pinion and gear speed, parallel 
| axis (rpm) 

Sac SAC allowable contact stress number 

Sat SAT allowable bending stress number 

SHP SHP shaft horsepower, output (hp) 

St SIGBE bending stress (psi) 
| SIG3P (E=epicyclic, P=parallel axis) 
! S. SIGHE contact stress (psi) 
3 SIGHP (E=epicyclic, P=paraliel axis) 
| T TORQE torgue (k in-lb) 
: TORQP (E=epicyclic, P=parallel axis) 

Wr WIE tangential tooth load (lb) 

WIP (E=epicyclic, P=parallel axis) 


II. COMMON BLOCK DETALLS 


The following provides information concerning the vari- 
ables in 2ach common block. The numbers in parentheses are 


she size of the array where applicable. 
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co G (FOR STRENGTH RATING) 


SFB : R®4 (2,2); service factor 


AKV 2 R&4; dynamic factor 
AKS 2 Rea; size factor 
AKM ¢ R24; load distribution factor 


AKO s R®GY (2) 3 overload factor 
SAT s R*4 (6)3 allowable bending stress number 


AKL 3 R®4 (2); life factor 


AKR 2 RXG (6); reliability factor 
AKT s RR4; temperature factor 
COMMON BLOCK AGM (FOR _ DOR R G) 


SFH : R®4 (2,2); service factor 


cv : R®4 (3); dynamic factor 

cs 2 RY; size factor 

CM 3 R®Y (2)3 load distribution factor 
CF 3s RAY; surface finish factor 

co 2 R®4 (2); overload factor 


SAC s R®&4 (6); allowable contact stress number 


cP s R*4; elastic properties factor 
CL s R*4 (2); life factor 

CH 2 R24; hardness factor 

cT s R&U3 temperature factor 

CR : R®G (6); reliability factor 
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COMMON BLOCK DESDAT (DESIGN PARAMETERS, INPUT) 
PWRIN : R™4 (2); source power input (hp) 
RPMIN : R™4 (2); source speed input (rpm) 
RPMOUT: R*®4; output shaft/propeller speed (rpm) 
DHELIX: R*& (3); helix angle (deg) 
HELIX =: R*4 (3)3 helix angle (rad) 
s R*& (3)3 transverse diametral pitch 
2 R&4 (3)3 normal diametral pitch 
R®4 (3)3 <Qansverse pressure angle (deg) 
RRQ (3)3 transverse pressure angle (rad) 
a™4 (3)3 normal pressure angle (dég) 
R*4 (3)3 Acrmal pressure angle (rad) 
number of different power sources 
arrangement code (1l=parallel axis, 
2=epicyclic) 
I*4 (3)3 epicyclic code (t=planetary, 2=Star) 
I*4 (3,2); hardness range code (1-6, see SUBR. 
AGMA) 
operational profile code (1=naval pro- 
file full power 5% max; 2=other, max 
power continuous) 
NPWRIN: I*4; number of power sources (inputs) 


ITPWRSR: power source code (t=turbine or motor, 


=multicylinder internai combustion 


engine) 
number of reduction stages 
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NPATH : I*4; number of power paths (1=single, 2=dual) 
NPLNT : I®4& (3)3 number of planet gears in epicyclic set 


NHELX 3; I*4; number of helicies (1=single, 2=double) 


COMMON BLOCK DESEPC (EPICYCLIC DESIGN PARAMETERS) 


4GOE : R®4; overall reduction ratio 


MGE : R®4 (3); stage reduction ratio 

RPMI =: R®Y (3); stage input speed (rpm) 

RPMPL : R*®4 (3); planet speed (rpm) 

ROMO 3 R®4 (3)3 stage output speed (rpm) 

OWRE =: R*4& (3);5 stage power split per planet (hp) 
(3)3 diameter of sun gear (in) 
(3); diameter of planet gears (in) 
(3) 3 diameter of ting gear (in) 
(3); facewidth (in) 
(3); durability geometry factor (sun/planet) 
(3); strength geometry factor (sun) 
(3)3 strength geometry factor (planet) 
(3)3 number of teeth, sun 
(3); number of teeth, planet 


(3)3 number of teeth, ring 


COMMON BLOCK DESPRL (PARALLE® AXIS DESIGN PARAMETERS) 
PWRFAC: R*4Y (2,3); stage power split factor 


MGOP <: R¥4 (2); overall reduction ratio 
MGP 3: R®Y (3,2); stage reduction ratio 
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DP 

DG 
FACED 
GEOMTI : 
GEOMJG: 
GEOMJP: 
NP 


NG 


TLPLE 3; 
UNTLDE: 
MPE 3 
KFCTRE: 
SIGHE ; 
SIGBE 

TORQE 

RPME : 
PDIAME: 


WGHTE 3; 


stage pinion and gear speed (rpm) 
stage power split per gear (hp) 
diameter of pinion (ina) 

diameter of gear (in) 

facewidth (in) 

duribility geometry factor 
strength geometry factor (gear) 
strength geometry factor (pinion) 
number of teeth, pinion 


number of teeth, gear 


(EPICYCLIC PARAMETERS, RESULTS) 
pitch line velocity (fpa) 
ft/d ratio (facewidth/sun diameter) 


center distance (theoretical) (in) 


tangential tooth load (1b) 


tooth load per in (lb/in) 
unit load (psi) 

mesh frequency (HZ) 

computed k-factor 

contact stress (psi) 

bending stress (psi) 

torgue (kK in=-1b) 

gear speeds (rpm) 

pitch diameters (in) 

gear set weight estimate (1b) 
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SPCWTE: R°4; specific weight (lb/hp) 
MTHE 3; I*®4 (3,3); tooth numbers 


ISIZES: I*4 (3); length, width, height estimates (in) 


COMMON BLOCK RESPRL (PARALLEL AXIS PARAMETERS, RESULTS) 
PLYP : R®4& (3,2); pitch line velocity (fpa) 


FBYD? ' (3,2); fd cratic (facewidth/pinion diameter) 
CDP (3,2); center distance (theoretical) (in) 
WTP ‘ (6,2); tangential tooth load (1b) 

TLPIP (6,2); tooth load per inch (l1b/in) 

UNTLDP (6,2); unit load (psi) 

MFP (3,2); mesh frequency (HZ) 

KFCTRP (6,2); computed k-factor 


SIGHP (3,2); contact stress (psi) 


SIGBP (6,2); bendiag stress (psi) 


TORQP (6,2); torque (k in-1lb) 

PDIAMP: R? (6,2); pitch diameters (in) 

SCDMIN: : minimum source center distance (in) 
SCDMAX: R*4; maximum source center distance (in) 
SHP : shaft horsepower, output (hp) 

aGHTP 3; gear set weight estimate (1b) 
SPCATP: specific weight (lb/hp) 

TRQOUT: R*4; torque, output (k in-lb) 

MTHP ;: I*4 (6,2); tooth numbers 


ISIZEP: I*4 (3); length, width, height estimates (in) 


APPENDIX ¢ 


REGAD SAMPLE RUNS 


This appendix contains samples of actual terminal ses- 
sions using REGAD. For the sake of brevity, only two con- 
plete sessions are included. However, a number of analysis 
and design runs were made using a full range of options and 
configurations, and they compared favorably to actual 
designs. The comparisons are not shown here due to the pro- 
prietary nature of the designs used for verification. The 


first example is an analysis run for a locked train, double 


reduction gear set with two different inputs. Following it, 


2s the results section from a design run using the identical 
parameters as the analysis run. The second example is a 
doubl2 reduction epicyclic gear set with the complete analy- 
sis session followed by the results section of a design run 
as before. The analysis and design sessions are identical 
with one exception. A seed for a random number generator is 
requested in the design option instead of diameters and 
facewidths as in the analysis option. For those cases where 
an infeasible design is generated, a message will alert the 
user and the program will continue. To obtain a feasible 


design, or just a different one, rerun the program and 


provide a different seed for the random number generator. 


This method was used on several occasions to obtain the 
desired results. Once a teasible design is obtained, the 
user can then use the analysis option to obtain a design 


that more closely suits his needs. 
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Figure 1: Single Reduction Parallel Axis Arrangements 
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Figure 2: Double Reduction Parallel Axis Arrangements 
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Figure 5: Triple Reduction Parallel Axis Arrangements 
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Figure 4: Single Reduction Epicyclic Arrangements 
(from Ref. 7) 


PRLDES PRLANL 


Program Description Option 


AGMA Constants Option 


PRLRES 


PRLSIZ 


PRLOUT 


EPCOOT 


Pigure 5: Flow Chart of the REGAD Package 
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Figure 6: 


External Tooth Dimensions 
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Figure ~ 


External Tooth Form Lavout 
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Figure 8: 


Fillet Center Location - 
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Figure 9: Fillet Center Location - Outside Base Circle 


Figure 10: Internal Tooth Dimensions 
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Figure 11: Internal Tooth Form Lavout 
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